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1 Abstract 

O ' 



The aim of this communication is to present in a concentrated form the main ideas of a 
method, developed by the author, for treating strongly nonequilibrium collective phenomena 
typical of the interaction of radiation with matter, as well as to give a survey of several ap- 
plications of the method. The latter is called the Scale Separation Approach since its basic 
techniques rely on the possibility of separating different space-time scales in nonequilibrium 
statistical systems. This approach is rather general and can be applied to diverse physical prob- 
^ ! lems, several of which are discussed here. These problems are: Superradiance of nuclear spins, 
' filamentation in resonant media, semiconfmement of neutral atoms, negative electric current, 

o 
o 
o 

£ 



and collective liberation of light. 



1 Introduction 

Strongly nonequilibrium processes that occur in statistical systems and involve their interac- 
tion with radiation are usually described by complicated nonlinear differential and integro- 
differential equations [1-3]. For treating these difficult problems, a novel approach has recently 
been developed [4-7] called the scale separation approach since its main idea is to formulate the 
evolution equations in such a form where it could be possible to separate several characteristic 
space-time scales. In many cases, different scales appear rather naturally being directly related 
to the physical properties of the considered system. 

The scale separation approach has been employed for solving several interesting physical 
problems related to strongly nonequilibrium processes occurring under the interaction of ra- 
diation with matter. As an illustration, the following phenomena are selected for this report: 
Superradiance of Nuclear Spins, Filamentation in Resonant Media, Semiconfinement of Neutral 
Atoms, Negative Electric Current, and Collective Liberation of Light. 

Since the scale separation approach makes the mathematical foundation for the following 
applications, its general scheme is described in Section 2. In Sections 3 to 7 concrete physical 
effects are briefly reviewed and the most important results are summarized. 
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2 Scale Separation Approach 



Because of the pivotal role of this approach for treating different physical problems, its general 
scheme will be presented here in an explicit way [4-7]. It is possible to separate the following 
main steps, or parts, of the approach. 

2.1 Stochastic quantization of short— range correlations 

When considering nonequilibrium processes in statistical systems, one needs to write evolution 
equations for some averages < Aj > of operators A^t) where t is time and i = l,2,...,N 
enumerates particles composing the considered system. For simplicity, a discrete index i is 
used, although everywhere below one could mean an operator A(fi, t) depending on a continuous 
space variable r*j. 

There is the well known problem in statistical mechanics consisting in the fact that writing 
an evolution equation for < Ai > one does not get a closed system of equations but a hierarchical 
chain of equations connecting correlation functions of higher orders. Thus, an equation for 

< Ai > contains the terms as J2j < A-iBj > with double correlators < A { Bj >, and the 
evolution equations for the latter involve the terms with tripple correlators, and so on. The 
simplest way for making the system of equations closed is the mean-field type decoupling 

< AiBj >^< Ai >< Bj >. When considering radiation processes, this decoupling is called 
the semiclassical approximation. Then the term £)■ < A^Bj > reduces to < A { > J2j < 
Bj >, so that one can say that < Ai > is subject to the action of the mean field Y,j < 
Bj >. The semiclassical approximation describes well coherent processes, when long-range 
correlations between atoms govern the evolution of the system, while short-range correlations, 
due to quantum fluctuations, are not important. However, the latter may become of great 
importance for some periods of time, for example, at the beginning of a nonequilibrium process 
when long-time correlations have had yet no time to develop. Then neglecting short-range 
correlations can lead to principally wrong results. 

To include the influence of short-range correlations, the semiclassical approximation can be 
modified as follows: 



where £ is a random variable describing local short-range correlations. It is natural to treat £ 
as a Gaussian stochastic variable with the stochastic averages 



where the second moment is defined so that to take into account incoherent local fluctuations. 
Since short-range correlations are often due to quantum fluctuations, the manner of taking them 
into account by introducing a stochastic variable £ can be called the stochastic quantization. 
Then the decoupling (1) may be termed the stochastic semiclassical approximation. This kind 




(1) 



« £ » = , « |£| 2 » = ]T I < Bj > 



(2) 
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of approximation has been used for taking into account quantum spontaneous emission of atoms 
in the problem of atomic superradiance. 



2.2 Separation of solutions onto fast and slow 

The usage of the stochastic semiclassical approximation makes it possible to write down a closed 
set of stochastic differential equations. The next step is to find such a change of variables which 
results in the possibility of separating the functional variables onto fast and slow, so that one 
comes to the set of equations having the form 

du ds 

— = f{e,u,s,Z,t) , — = e g{e,u,s,£,t) , (3) 

where e <C 1 is a small parameter, such that 

lim f ^ , lim e q = . (4) 

As is evident, dealing with only two functions, u and s, and one small parameter e is done just 
for simplicity. All procedure is straightforwardly applicable to the case of many functions and 
several small parameters. 

From Eqs. (3) and (4) it follows that 

lim ^ , lim ^ = , (5) 

which permits one to classify the solution u as fast, compared to the slow solution s. In turn, 
the slow solution s is a quasi-invariant with respect to the fast solution u. 

The above classification of solutions onto fast and slow concerns time variations. In the 
case of partial differential equations, one has, in addition to time, a space variable r. Then the 
notion of fast or slow functions can be generalized as follows [8,9]. Let feV, with mesV = V, 
and t G [0, T], where T can be infinite. Assume that 

* f du I -> 

lin ^ < T7 / ^7 rff> ^ ' Vudt^>^0, (6) 

e^O V Jv (Jt e-+0 1 Jo 

while 

lim < — / — dr > = , lim < — / Vs dt > = . (7) 



e^O V JV dt ~ e^O T JO 

Then the solution u is called fast on average, with respect to both space and time, as compared 
to s that is slow on average. In such a case s is again a quasi-invariant as compared to u. In 
general, it may, of course, happen that one solution is fast with respect to time but slow in 
space, or vice versa, when compared to another function. The notion of quasi-invariants with 
respect to time is known in the Hamiltonian mechanics where they are also called adiabatic 
invariants. Here this notion is generalized to the case of both space and time variables [8,9]. 
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2.3 Averaging method for multifrequency systems 



After classifying in Eqs. (4) the function u as fast and s as slow, one can resort to the Krylov- 
Bogolubov averaging technique [10] extended to the case of multifrequency systems. This is 
done as follows. 

Since the slow variable s is a quasi-invariant for the fast variable u, one considers the 
equation for the fast function u, with the slow one kept fixed, 

8X 

— = f(e,X,z,U)- (8) 

Here z is treated as a fixed parameter. The solution to Eq. (8), that is 

X — X(e, z, £, t) , z = const , (9) 

has to be substituted into the right-hand side of the equation for the slow function, and for 
this right-hand side one defines the average 

g(e, z) = < - T g(e, X(e, z, £, t), z, £, t) dt > , (10) 
r Jo 

in which r is the characteristic oscillation time of the fast function. In many cases, it is possible 
to take t — > oo, especially when the period of fast oscillations is not well defined [2]. Then one 
comes to the equation 

dz _ 

— = eg(e,z) (11) 

defining a solution 

z = z(e,t). (12) 

Substituting the latter into X, one gets 

y(e,Z,t)=X(e,z(e,t),Z,t). (13) 

The pair of solutions (9) are called the generating solutions since these are the first crude 
approximations one starts with. More elaborate solutions are given by Eqs. (12) and (13) 
which are termed guiding centers. 

Notice two points that difference the case considered from the usual averaging techniques. 
The first point is that in Eq. (8) the small parameter e is not set zero. And the second difference 
is in the occurrence of the stochastic average in Eq. (10). Leaving e in Eq. (8) makes it possible 
to correctly take into account attenuation effects, as will be shown in applications. 



2.4 Generalized expansion about guiding centers 

Higher-order corrections to solutions may be obtained by presenting the latter as asymptotic 
expansions about the guiding centers (12) and (13). To this end, /c-order approximations are 
written as 

k 

u k = y(e, f , f) + 2 Vn(e, f , t) e n , 

n=l 
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k 

s k = z(e,t) + Y,Zn(e,t,t)e n . (14) 

n=l 

Such series are called generalized asymptotic expansions [11] since the expansion coefficients 
depend themselves on the parameter e. The right-hand sides of Eqs. (3) are to be expanded 
similarly to Eq. (14) yielding 

k 

f(e, u k , s fc , f , t) ~ /(e, y, z, £, i) + ^ /n(e, f , £™ (15) 

n=l 

and an equivalent expansion for g. These expansions are to be substituted into Eqs. (3) with 
equating the like terms with respect to the powers of e. In the first order, this gives 

= A( e »f>*) -9(e,z)Xi(e,Z,t) , < ^- = g{e,y,z,$;,t)-g{e,z), (16) 

where 

d 

Xi(e, f , t) = — X(e, z,£,t) , z = z(e, t) . 
For the approximations of order n > 2, one gets 

^f = M £ >t> t )> ^ = 9n{e^t). (17) 

The functions f n and g n depend on yi,y2, ■ ■ ■ ,y n and on zi,z 2 , . . . , z n (see for details [6,7]). 
But it is important that the dependence on y n and z n is linear. Therefore all equations (16) 
and (17) are linear and can be easily integrated. Thus, the approximants (14) are defined. 
Each /c-order approximation can also be improved by invoking the self-similar summation of 
asymptotic series [12-18]. 



2.5 Selection of scales for space structures 

The solutions of differential or integro-differential equations in partial derivatives are often 
nonuniform in space exhibiting the formation of different spatial structures. Also, it often hap- 
pens that a given set of equations possesses several solutions corresponding to different spatial 
patterns or to different scales of such patterns [3]. When one has a set of solutions describing 
different possible patterns, the question arises which of these solutions, and respectively pat- 
terns, to prefer? This problem of pattern selection is a general and very important problem 
constantly arising in considering spatial structures. In some cases this problem can be solved 
as follows. 

Assume that the obtained solutions describe spatial structures that can be parametrized 
by a multiparameter b, so that the /c-order approximations Uk(b,t) and Sk(b,t) include the 
dependence on b whose value is however yet undefined. To define b, and respectively the 
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related pattern, one may proceed in the spirit of the self-similar approximation theory [12- 
14], by treating b as a control function, or a set of control functions if b is a multiparameter. 
According to the theory [12-14], control functions are to be defined from fixed-point conditions 
for an approximation cascade, which is to be constructed for an observed quantity. For the 
latter, one may take the energy which is a functional E[u, s] of the solutions. In experiments, 
one usually measures an average energy whose /c-order approximation writes 

£*(&) = «- [ T E[u k (b,t),s k (b,t)]dt^> , (18) 

T JO 

where r is a period of fast oscillations. For the sequence of approximations, {E k (b)}, it is 
possible to construct an approximation cascade [12-14] and to show that its fixed point can be 
given by the condition 

^( b ) = , (19) 

from which one gets the control function b = b k defining the corresponding pattern. According 
to optimal control theory, control functions are defined so that to minimize a cost functional. 
In this case, it is natural to take for the latter the average energy (18). Therefore, if the fixed- 
point equation (19) has several solutions, one may select of them that one which minimizes the 
cost functional (18), 

E k {b k ) =absmin£ fc (&) . (20) 

b 

Equations (19) and (20) have a simple physical interpretation as the minimum conditions for 
the average energy (18). However, one should keep in mind that there is no in general such a 
principle of minimal energy for nonequilibrium systems [3]. Therefore the usage of the ideas from 
the self-similar approximation theory [12-14] provides a justification for employing conditions 
(19) and (20) for nonequilibrium processes. 

In the following sections a brief survey is given of several physical examples the scale sepa- 
ration approach has been applied to, and the main results are formulated. 



3 Superradiance of Nuclear Spins 

A system of neutral spins in an external magnetic field , prepared in a strongly nonequilib- 
rium state and coupled with a resonance electric circuit, displays rather nontrivial relaxation 
behaviour somewhat similar to that of an inverted system of atoms. This is why the optical 
terminology, such as superradiance, has been used for describing collective relaxation processes 
in nonequilibrium nuclear magnets [5,6,19,20]. 

For a system of nuclear spins interacting through dipole forces the evolution equations can 
be derived [5,6] for the averages 

i N i N 

^^£<*r>, 8 = -52<s;>, (21) 
jv i=\ jv i=i 
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in which N is the number of spins, angle brackets mean statistical averaging, S~ is a lowering 
spin operator, and Sf is the ^-component of a spin operator. Following the ideology of the 
scale separation approach, local fluctuating fields are presented by stochastic variables £ an d 
£. In this way, one comes to the evolution equations for the transverse spin variable 

du 

— = i{uj - ^0 + ^72)^-^(73^ + s ( 22 ) 

and the longitudinal average spin 

d s 1 1 

Tt = ^(lsh + On* - -(l 3 h + C)u - 7i(* - ■ (23) 



It is also convenient to consider the equation 

d_ 
di 



u\ 2 = -272M 2 - i(7 3 /i + £)su* + i(7 3 /i + C)su . (24) 



In equations (22)-(24) dimensionless units are used for the resonator magnetic field h satisfying 
the Kirchhoff equation 

dh r* d 

— + 2 l2 h + u 2 / h(t') di! = -2a - (u* + u)+ 73/ • (25) 
at Jo at 

Here u>o is the Zeeman frequency of spins in an external uniform magnetic field, u is the 
resonator natural frequency, 71 and 72 are the spin-lattice and spin-spin relaxation parameters, 
respectively, 73 is the resonator ringing width, ( is a stationary spin polarization, a is the 
coupling between spins and the resonator, and / is an electromotive force. The random local 
fields are defined as Gaussian stochastic variables with the stochastic averages 

« e » = « Kl 2 » = 7 2 * , (26) 

where 72 is the inhomogeneous dipole broadening. 

There are the following small parameters in the system: 

— < 1 , — < 1 , — < 1 , — < 1 , 

LU UJq UJq UJ 

— < 1 , (A = uj - uj ) . (27) 

UJq 

This makes it admissible to classify the functions u and h as fast, while s and \u\ 2 as slow, and 
to apply the method of Section I. The behaviour of solutions to Eqs. (22)-(25) depends on 
initial conditions for u(0), and s(0), on the existence of an electromotive driving force f(t), on 
the pumping related to the parameter (, and on the value of the effective coupling parameter 

2 P /4 UJQ 

g = n 77 , 28 

n 72 uj 
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in which r] is a filling factor; p, spin density; and [i n is a nuclear magnetic moment. 

The first interesting result is that the electromotive force does not influence much macro- 
scopic samples [5,6] since the corresponding correlation time is proportional to N, that is, the 
effective, interaction strength of an electromotive force with the spin system is proportional 
to N^ 1 . This shows, in particular, that the role of the thermal Nyquist noise for starting the 
relaxation process is negligible. The main cause triggering the motion of spins leading to coher- 
ent self-organization is the presence of nonsecular dipole interactions [5,6,19]. The latter result 
gives an answer to the problem, posed by Bloembergen and Pound [21]: What is the origin of 
self-organized coherent relaxation in spin systems? 

All possible regimes of nonlinear spin dynamics have been analysed [5,6,19,20]. When the 
nonresonant external pumping is absent, that is ( > 0, there are seven qualitatively differ- 
ent transient relaxation regimes: free induction, collective induction, free relaxation, collective 
relaxation, weak superradiance, pure superradiance, and triggered superradiance [6] . In the pres- 
ence of pumping, realized e.g. by means of dynamical nuclear polarization directing nuclear 
spins against an external constant magnetic field, one has ( < 0. Then three dynamical regimes 
can be observed, depending on the value of ( with respect to the pumping thresholds 

Ci = "-, C2 = --(l + ^) , (29) 
9 9 \ 2 72y / 

where 7J 1 is an effective pumping rate. 

Two stationary points can exist for the slow solutions s and w, where 

w = \u\ — 2 — s . 



LU 



These fixed points are 



si = C , w\ = , 

1 * 7i(l + Q ( or,\ 

S 2 = , W 2 = 5 • ( 30 ) 

9 9 2 l2 

When d < ( < 0, then the first fixed point is a stable node and the second one is a saddle 
point. For ( = (1, both stationary points merge together, being neutrally stable. After the 
bifurcation at the value ( = (1, in the region ( 2 ^ C < Ci ? the first fixed point looses its stability 
becoming a saddle while the second fixed point becomes a stable node. Finally, when ( < ( 2 , 
the second fixed point transforms to a stable focus, and the first one is, as earlier, a saddle 
point. 

In this way, there are three qualitatively different lasting relaxation regimes induced by the 
pumping. The first one is a monotonic relaxation to the first stationary solution with practically 
no coherence, w{ = 0. The second regime is a monotonic relaxation to the second stationary 
solution with a nonzero coherence, 7^ 0. And the third regime is that of pulsing relaxation 
to the coherent stationary point. Note that the pumping rate 7^ can be larger than 72, so that 
u>2 can reach the order of unity. The three lasting relaxation regimes occurring in the presence 
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of pumping can be called, respectively, incoherent monotone attenuation, coherent monotone 
relaxation, and coherent pulsing relaxation. 



4 Filamentation in Resonant Media 

In optical resonant media there appear space structures when the radiation wavelength is much 
less than the characteristic sizes of the laser system [22]. There are two principally different 
types of spatial structures in laser media, one corresponding to low Fresnel numbers [23-28] 
and another type corresponding to high Fresnel numbers [29-35], with a transition occurring 
around F ~ 10. Similar effects are observed in photorefractive media [36-38]. Such structures 
are described by nonlinear differential equations in partial derivatives. The general problem in 
dealing with these equations is the nonuniqueness of their solutions each of which corresponds 
to a particular spatial structure [3]. The related problem of pattern selection can be treated by 
the method of subsection 2.5. Here this is illustrated by the theory of filamentation in optical 
resonant media [39-42] at high Fresnel numbers. 

The Hamiltonian for a system of resonant atoms can be written as 

h = \y: + °D -\z (? • ^Vr + °td- 4i) - 

i i 

-lY,(d*-E±a-+atd.E^ , (31) 

where the standard notation is used for the Pauli matrices of and erf; d is a transition dipole 
for a transition with the frequency ojq, the electric fields 



E~ = Eoe^-^ (k = ^ 

E7 j = !$- fUj x (dx n %j ) e ikor -a- (k = ^) (32) 
correspond to the laser mode and to the reradiated field, respectively, and 

I -» I -, f'ij -, -, -, 

'I j — I ' IJ I ) 'Hj — 5 ' ij — ' I 1 J ■ 

r ij 

The resonant medium has cylindrical shape of radius R and length L. The transition wavelength 
A is such that 

A R , „ 

- « 1 , - « 1 . (33) 

It is convenient to pass to a continuous space variable f by transforming the sums in integrals 
according to the rule 

22fi = PJ f( f ) df \P=y 
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Then the evolution equations for the statistical averages 



u(r,t) =< a (f,t) > , s(f,t) =< a z (f,t) > (34) 
satisfy partial integro-differential equations. Because of the inequalities 

To I A I 

— < 1 , — < 1 , (35) 

where A = uj — ojq is the detuning parameter, the function u is fast in time as compared to the 
slow function s. 

The solutions to the evolution equations describe a bunch of filaments aligned along the z 
axis which is the axis of the sample. Each filament has a radius Rf and is centered at a point 
{x n ,y n } in the transverse cross-section. The location of the filament centers is distributed 
chaotically. Thus, the solutions may be presented as expansions over filaments in the form 

Nf / \ 

u(r,t) = ]T u n (t)e ikz 6 (R f - y/(x - x n ) 2 + (y - y n f ) , 

n=l V ' 



N f , 

s (r,t) = J2 s n(t) e (R f - ^( x - Xn y + ( y 

n=l V 



y n ) 2 ) , (36) 



where Nf is the number of filaments and ©(•) is a unit-step function. The number of filaments 
is related to the filament radius Rf and the pumping characteristique 

C{t) = \ J s(f,t) df , (37) 

which yields 

^ = i( 1 + «(|) 2 - < 38 > 

The filament radius Rf can be defined according to the procedure of subsection 2.5. To this 
end, Rf is considered as a control function parametrizing the filamentary space structure. 
It is possible to construct the average energy (18) corresponding to the Hamiltonian (31). 
Minimizing this energy functional with respect to Rf gives 



R f = 0.22 VU . (39) 
Then the number of filaments (38) can be written as 

N f = 3.3 (1 + C)F ( F = ! Xr)> ( 4 °) 

where F is a Fresnel number. 

The predictions of the theory [39-42] have been found to be in a good agreement with 
measurements, as has been confirmed in a series of experiments [29-34] with different lasers. 



10 



5 Semiconfinement of Neutral Atoms 



Dynamics of neutral atoms in nonuniform magnetic fields concerns problems of current exper- 
imental and theoretical interest, especially with regard to atoms in magnetic traps, where the 
atoms can be cooled down to experience the Bose-Einstein condensation [43,44]. The motion 
of confined atoms is usually described by means of the adiabatic approximation assuming that 
the atom spins are permanently aligned along the local magnetic field and, thus, adiabatically 
follow its direction. To study the more general case, when atoms are permitted to escape from 
a trap, one has to invoke a more refined approximation, such as the scale separation approach 
described in Section 2. 

The motion of neutral atoms in magnetic fields can be presented by the semiclassical equa- 
tions for the quantum-mechanical average of the real-space coordinate, R = {R a }, and for the 
average S = {S a } of the spin operator, with a — x,y, z. The first equation writes 

d 2 R a Ho x dB 

-^=m S -dR a +li ^ (41) 

— * 

where /x is magnetic moment, m is mass of an atom, B is a magnetic field, and 7£ a is a collision 
term. The equation for the average spin is 

f^SxB. (42) 

The total magnetic field B = B\ + B 2 consists of two terms. One is the quadrupole field 

B x = B[ (R x e x + R y e y + XR z e z ) , (43) 
in which A is the anisotropy parameter. The second term is a transverse field 

— # 

B 2 = B 2 (e x cos cut + e y sin tot) . (44) 
The characteristic frequencies 



Z^o-BiV^ 2 _ HqB : 



- s (SJ • (45) 

where Ro = B 2 /B[, satisfy the inequalities 

^ « 1 , - « 1 . (46) 

UJ 2 UJ 2 

Because of the latter, the spin variable S has to be classified as fast, compared to the slow 
atomic variable R. 
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The collision term in Eq. (41) contains a collision rate 7 and a random collision variable 
£ a (t) defined by the stochastic averages 

« Ut) > = , < &(*)&(0 > = ZD a 5 aP 5{t - t') , (47) 

where D a is a diffusion rate. 

The semiconfining regime of motion can be realized by preparing for the spin variable 
nonadiabatic initial conditions 

S° x = S° y = , S° z = S ? , (48) 

which can be done e.g. by means of an external pulse at t = 0. Then it is possible to show 
[45-47] that the motion of atoms becomes axially restricted by the value 



z m Ro 

such that 



mm t R z (t) (XS>0) 
max t R z (t) (XS<0), 

3z 2 

z ™ = z » ~~ 2>dui ■ (49) 

Atomic collisions do not disturb the semiconfined motion provided that temperature T is suf- 
ficiently low satisfying the condition 

« ! , (50) 



1 



where p is the density of particles and a is a scattering length. 

The semiconfining regime of motion makes it possible to form well-collimated beams of 
neutral particles by means of only magnetic fields. This regime can be employed for creating 
coherent beams of Bose atoms from atom lasers. 



6 Negative Electric Current 

Semiconductors with nonuniform distribution of charge carriers can exhibit rather unusual 
transport properties. For example, in a sample biased with an external constant voltage, the 
transient effect of negative electric current can happen [8,9]. 

Transport properties of semiconductors are usually described by the semiclassical drift- 
diffusion equations. In what follows a plane device is considered and all quantities are expressed 
in dimensionless form, so that the space variable is x G [0, 1]. The continuity equation writes 
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where p(x,t) is a charge density; E(x,t) is the electric current; Di, and Tj are mobility, 
diffusion coefficient, and relaxation time, respectively, for the carriers of type i. The Poisson 
equation is 

f£ = «*!>. (52) 

At the initial time, the distribution of charge carriers is nonuniform, given by 

Pi (x,0) = f i (x). (53) 
The sample is biased with an external constant voltage, which means that 

f 1 E(x, t) dx = 1 . (54) 
Jo 

The total electric current through the semiconductor sample is 

J(t) = /^(M) dx, (55) 

where the density of current 

Owing to the voltage integral (54), one has 

f 1 4- E(x,t) dx = . (56) 
Jo at 

It is also possible to show that 

lim - [ T ^-E(x, t) dt = . (57) 

t^oo T Jo OX 

This means that the function E can be considered as slow on average in time and space. Then, 
treating E as a quasi-invariant, one may find the solutions to Eqs. (51) and (52) in order to 
analyse their general space-time behaviour. 

Negative electric current can appear only when the initial charge distribution is essentially 
nonuniform. If this initial distribution forms a narrow layer located at the point x — a, then 
the current (55) becomes negative for some short time close to t — 0, if one of the following 
conditions holds true: 

1 1 



a < - 



2 4nQ 



(«>*)• 



a > I + t-U (q < ) . ( 58 ) 



2 4vr|Q| V 2tt 
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where 



Q = J2Qi ^ Qi = Pi(x, 0) dx 



(59) 



The effect of the negative electric current can be employed for various purposes, as is discussed 
in Refs. [8,9]. For instance, when the initial charge layer is formed by an ion beam irradiating 
the semiconductor sample, then the location a corresponds to the mean free path of these ions. 
In this case, by measuring the negative current J(0), one can define this mean free path 



1 

a= 2 



4ttQ 



AO) 



(60) 



This formula is valid for both positive and negative values of Q. 



7 Collective Liberation of Light 

One more interesting physical effect that has been described using the scale separation approach 
is collective liberation of light [48]. Consider an ensemble of resonant atoms which are doped 
into a medium with well developed polariton effect [49] , when in the spectrum of polariton states 
there is a band gap. If an atom with a resonance frequency inside the polariton gap is placed into 
such a medium, the atomic spontaneous emission is suppressed, which is called the localization 
of light [50] . However, a system of resonant atoms inside the polariton gap can radiate if their 
coherent interaction is sufficiently strong. Thus the suppression of spontaneous emission for a 
single atom can be overcome by a collective of atoms radiating coherently. Conditions when 
such a collective liberation of light can arise and the dynamics of this liberation are analysed 
inRef. [48]. 

In conclusion, a general method has been developed for treating strongly nonequilibrium 
processes in statistical systems. This method, called the scale separation approach, is especially 
useful for describing collective phenomena in the interaction of radiation with matter. To 
emphasize the generality of the approach, it is illustrated here by several different physical 
examples. The common feature of all considered systems is that their evolution is described 
by nonlinear differential or integro-differential equations. Such equations, as is known, are 
difficult to solve. The scale separation approach makes it possible to find accurate approximate 
solutions. The accuracy of the latter has been confirmed by numerical calculations and by 
comparison with experiment, when available. Using this approach several interesting physical 
problems have been solved and new effects are predicted. 
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